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Abstract — In this paper, a general criterion for space time block 
codes (STBC) to achieve full-diversity with a linear receiver is 
proposed for a wireless communication system having multiple 
transmitter and single receiver antennas (MISO). Particularly, 
the STBC with Toeplltz structure satisfies this criterion and 
therefore, enables full-diversity. Further examination of this 
Toeplitz STBC reveals the following Important properties: a) 
The symbol transmission rate can be made to approach unity, 
b) Applying the Toeplltz code to any signalling scheme having 
nonzero distance between the nearest constellation points results 
in a non-vanlshlng determinant. In addition. If QAM is used as 
the signalling scheme, then for Independent MISO flat fading 
channels, the Toeplitz codes Is proved to approach the optimal 
diverslty-vs-multiplexing tradeoff with a ZF receiver when the 
number of channel uses Is large. This is, so far, the first non- 
orthogonal STBC shown to achieve the optimal tradeoff for such 
a receiver. On the other hand, when ML detection is employed 
in a MISO system, the Toeplltz STBC achieves the maximum 
coding gain for Independent channels. When the channel fading 
coefficients are correlated, the inherent transmission matrix in 
the Toeplltz STBC can be designed to minimize the average worst 
case pair-wise error probability. 

Index Terms — Full Diversity, Linear Receiver, MISO, ML 
detection. Non-vanishing determinant. Optimal dlversity-vs- 
multiplexing tradeoff, STBC, Toeplltz 



I. Introduction 

THE recent arrival of the Information Age has created 
an explosive demand for knowledge and information 
exchange in our society. This demand has triggered off an 
enormous expansion in wireless communications in which 
severe technical challenges, including the need of transmitting 
speech, data and video at high rates in an environment rich 
of scattering, have been encountered. A recent development 
in wireless communication systems is the multi-input multi- 
output (MIMO) wireless link which, due to its potential in 
meeting these challenges caused by fading channels together 
with power and bandwidth limitations, has become a very 
important area of research. The importance of MIMO com- 
munications lies in the fact that they are able to provide a 
significant increase in capacity over single-input single-output 
(SISO) channels. Existing MIMO designs employ multiple 
transmitter antennas and multiple receiver antennas to exploit 
the high symbol rate provided by the capacity available in 
the MIMO channels. Full symbol rate is achieved when, on 
average, one symbol is transmitted by each of the multiple 
transmitter antennas per time slot (often called a "channel 
use")- In the case of M transmitter antennas, we will have 
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an average of M symbols per channel use (pcu) at full 
rate. Furthermore, to combat fading and cross-talk, MIMO 
systems provide different replicas of transmitted symbols to 
the receiver by using multiple receiver antennas with sufficient 
separation between each so that the fading for the receivers 
are independent of each other. Such diversity can also be 
achieved at the transmitter by spacing the transmitter antennas 
sufficiently and introducing a code for the transmitted symbols 
distributed over transmitter antennas (space) and symbol pe- 
riods (time), i.e., space-time coding [l]-[4]. Full diversity is 
achieved when the total degree of freedom available in the 
multi-antenna system is utilized. 

Over the past several years, various space-time coding 
schemes have been developed to take advantage of the MIMO 
communication channel. Using a linear processor, orthogonal 
space-time block codes [2], [3], [5]-[8] can provide maximum 
diversity achievable by a maximum likelihood detector. How- 
ever, they have a limited transmission rate [8]-[ll] and thus, 
do not achieve full MIMO channel capacity [12]. Linear dis- 
persion codes have been proposed in [13] for which each trans- 
mitted codeword is a linear combination of certain weighted 
matrices maximizing the ergodic capacity of the system. 
Unfortunately, good error probability performance for these 
codes is not strictly guaranteed. To bridge the gap between 
multiplexing and diversity, a linear dispersion code design has 
been proposed using frame theory [14] that typically performs 
well both in terms of ergodic capacity and error performance, 
but full diversity still cannot be guaranteed. Thus far, with 
the exception of the orthogonal STBC, all existing STBC 
are designed such that full diversity can only be achieved 
when the ML detector is employed. Recent research [15]- 
[19] based on number theory has shown that employing a ML 
receiver, it is possible to design linear space-time block codes 
and dispersion codes which are full rate and full diversity 
without information loss. The major concern on these designs 
is that the coding gain vanishes rapidly as the constellation size 
increases. Therefore, designs of full-rate, full-diversity space- 
time codes with non-vanishing coding gain have drawn much 
attention [20]-[32] since such structured space-time codes 
could achieve the optimal diversity-vs-multiplexing tradeoff 
developed by Zheng and Tse [33]. However, most available 
STBC possessing these properties are for ML receivers only. 

In this paper, we consider a coherent communication system 
equipped with multiple transmitter antennas and a single 
receiver antenna, i.e., a MISO system. These systems are often 
employed in mobile communications for which the mobile 
receiver may not be able to support multiple antennas. The 
highest transmission rate for a MISO system is unity, i.e.. 
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one symbol pcu. For such a MISO system with ML receivers, 
rate-1 and full diversity STBC have been proposed by various 
authors [29], [34]-[36]. In this paper, however, we consider 
such a MISO system equipped with linear receivers for which 
we propose a general criterion for the design of a full-diversity 
STBC. In particular, we introduce the Toeplitz STBC as a 
member of the family of the full diversity STBC. It should be 
noted that the Toeplitz structure has already been successfully 
employed as a special case of the delay diversity code (DDC) 
[37]-[40] applied to MIMO systems having outer channel cod- 
ing and ML detection. Here, we extend its application to the 
construction of STBC in a MISO system by having a Toeplitz 
coding matrix cascaded with a beamforming matrix. We show 
that the Toeplitz STBC has several important properties which 
enable the code, when applied to a MISO system with a 
linear receiver, to asymptotically achieve unit symbol rate, to 
possess non-vanishing determinants for signal constellations 
having non-zero distance between nearest neighbours, and to 
achieve full diversity [41] accomplishing the optimal tradeoff 
of diversity and multiplexing gains [33]. 

On the other hand, we also consider the MISO system 
in which the channel has zero mean and fixed covariance 
known to the transmitter For such MISO systems, sacrific- 
ing the transmission rate by repeating the transmitted sym- 
bols, and employing maximum ratio combining together with 
orthogonal space-time coding, an optimal precoder can be 
designed [42], [43] by minimizing the upper bound of the 
average symbol error probability (SEP). Here in this paper, we 
apply the Toeplitz STBC to such a MISO system. Maintaining 
rate one and full diversity, we present a design that minimizes 
the exact worst case average pair-wise error probabiUty when 
the ML detector is employed at the receiver 

II. MISO System Model and Properties of the 
Channel Matrix 

Consider a MIMO communication system having M trans- 
mitter antennas and Mji receiver antennas transmitting the 
symbols {si}, £ = 1, . . . ,L which are selected from a given 
constellation, i.e., se G S. To facilitate the transmission of 
these L symbols through the M antennas in the N time 
slots (channel use), each symbol is processed by an A^ x M 
coding matrix A^, and then summed together, resulting in 
an N X M STBC matrix given by X = J2e=i ^i-^i where 
the (nm)th element of X represents the coded symbol to be 
transmitted from the mth antenna at the 7ith time slot. These 
coded symbols are then transmitted to the receiver antennas 
through flat-fading path coefficients which form the elements 
of the M X Mji channel matrix H. The received space-time 
signal, denoted by the N x Mji matrix Y, can be written as 



Y = XH + H 



(1) 



where H is the N x Mr additive white space-time noise matrix 
whose elements are of complex circular Gaussian distribution 
CAA(0,1). 

Let us now turn our attention to a MISO wireless commu- 
nication system which is a special case of the MIMO system 
having M transmitter antennas and a single receiver antenna. 



Just as in the MIMO system, the transmitted symbols sg, 
i = 1, • • • , i in the MISO system are coded by Hnear N x Ad 
STBC matrices Ai which are then summed together so that 



X = ^A, 



(2) 



where L is the total number of symbols to be transmitted If 
L ~ N, the system is at full-rate (rate-one). At the time slot 
77,, the nth row of the coding matrix X feeds the M coded 
symbols to the M antennas for transmission. Each of these 
transmitter antennas is linked to the receiver antenna through 
a channel path coefficient hm, ?n = 1, • • • , M. At the receiver 
of such a system, for every N time slots (77, = 1, • • • , N), we 
receive an iV-dimensional signal vector y = [j/i 7/2 • • • Vn]'^ 
which, as a special case of Eq. ([T]), can then be written as 



y = xh + | 



(3) 



where h = [hi,- ■ ■ , hnj]'^ is an il/x 1 channel vector assumed 
to be circularly symmetric complex Gaussian distributed with 
zero-mean and covariance matrix S, and ^ is an iV x 1 noise 
vector assumed to be circularly symmetric complex Gaussian 
with covariance (t'^In- Putting Eq. ^ into Eq. (O, writing the 
symbols to be transmitted as a vector and aligning the code- 
channel products to form the new channel matrix we can write 



n = (Aih Ash 



Aih) and s = [si S2 • • • sl] (4) 



the received signal vector can now be written as 

y = ^s + ^ 



(5) 



In this paper, we emphasize on the application of linear 
receivers for the MISO system in Eq. (|5]). In the following, 
we will derive a condition on the equivalent channel Ti. that 
renders full-diversity when the signals are received by a linear 
receiver First, we present the following properties of the 
equivalent channel matrix Ti,: 

Property 1: Suppose the equivalent channel Ti. in Eq. (HI 
is such that 7t 7i is non-singular for any nonzero h. Then 
we have the following inequality: 



C„ 



\2L 



< det 



n"n 



<c„ 



\2L 



(6) 



where Cmin and C,, 
h. 



are positive constants independent of 



Proof: Since h is nonzero, we normalize the L x L 
matrix 7i 7i hy dividing each of its elements with ||h||^, 
i.e., Ti. Ti. ~ ||h|p]HI, where H is the normalized matrix with 
the ijih element being equal to 

[H]y- = ^i^Af A, ^ 7, j = 1, 2, • • • , L 



The determinant of positive semi-definite (PSD) matrix H is 
continuous in a closed bounded feasible set {h : ||hjp = 1} 
where h = ttt- It has the maximum and minimum values 

l|h| 

that are denoted by C„iax and Cmin respectively. Now, since 
7i 7^ is non-singular for any nonzero h, its determinant is 
^mm < Cinax and Eq. ^ holds. D 



positive. Therefore, < C^ 
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The following example serves to illustrate the above prop- 
erty. 

Example 1: Consider the following channel matrix H 
^hi 0' 

The determinant of matrix Ti, Ti, can be written 




Aeti-H^H) 









Since 

■U4 = sin 6, and Eq. (|7]i becomes 



1, we can define '—^ 



\hi? \h2? 

I|h||2||h!|2 
\hl\ _ 



dciiU^n) 









sm 
1 



SIT? {29) 



(7) 
cos 6', and 

(8) 



I- \ SIT? {20) is continu- 



It is obvious that the function f{6) — ^ ^ . 
ous in a closed bounded set. The minimum and maximum of it 
can be easily obtained as Cmin = |; C^ax = 1- Both values 
are constants and are independent of the random channel. 
Thus, the determinant of the channel matrix is bounded 
by |||h||4 < dct (t^^-h) < ||h!|4. D 

Property 2: If Ti. Ti. is non-singular for any nonzero h, 
then the diagonal elements of [Ti Ti.]^^ satisfies the follow- 
ing inequality 

-1 



(n^'n) 



>CoM 



(9) 



for £ = 1, 2, • • • , L where Co is a constant independent of h.H 
Proof: From the matrix inversion algorithm [44], we have 



(n"n 



dct 



n"n 



det 



I Hg Hi 



(10) 



where T-ii is the matrix obtained by deleting the £th column 
vector, Afh, from Ti. We notice that the matrix T-ii Ti-i is 
still PSD and therefore satisfies the right side inequality of 
Eq. (|6]l having an upper bound denoted by CfmaxUhH^'^^^^. 
Applying the lower bound of Eq. (|6]l to the numerator and the 
upper bound to the denominator of Eq. (fTOl i. we obtain 



{n"H) 



> 



Cn 



\2L 



U 






L-min/Ccj 



with Qmax 



h||2(i-l) 

h||'>Coi|h||2 



(11) 



max{C| 



£ max 1 ^ 



where Co = 

1,2,--.,L}. D 

Properties [Hand ID are of fundamental importance to the design 
of full diversity STBC for a MISO system employing a linear 
detector This will be presented in the following section. 

III. Diversity Gain of STBC for a MISO System 
Employing a Linear Receiver 

Let us first review the concept of diversity gain with 
reference to a MIMO system. Consider the MIMO system in 



Eq. dU equipped with a maximum likelihood (ML) detector. It 
is well-known that an upper bound for the average pair-wise 
error probability is given by [45] 



P(s 



< 



< 



1 



dct I 



^M 



Mn^" 



8M 

-Mr 



X^(e)X(e) 



-Mn 



\8MJ 



-rMa 



(12) 



where p = E tr[X^X]/E tr[H^H] is the SNR, e = s - s' 
with s, s' G 5 is the error vector, r {< M) is the rank, 
and {Am},m ~ !,■■■ ,r are the non-zero eigenvalues of 
the matrix X^(e)X(e). The middle part of Eq. (fTZt is the 
Chernoff bound, which at high SNR, can be further tightly 
bounded by the right side. For a given Mji, two factors dictate 
the minimization of this bound on the right side of Eq. (fT2b : 

a) The Rank of X^(e)X(e): The exponent rMu of the 
second term governs the behaviour of the upper bound 
with respect to SNR and is known as the diversity gain. 
To keep the upper bound as low as possible, we should 
make the diversity gain as large as possible. Full diversity 
is achieved when r = M, i.e., X(e) is of full column 
rank. This implies that the diversity gain achieved by an 
ML detector depends on e, which is decided by the type 
of signalling. 

b) The Determinant of X^(e)X(e): The first term consists 
of the product of the non-zero eigenvalues of X^(e)X(e) 
and is called the coding gain. For X^(e)X(e) being full 
rank, this product is its determinant the minimum value of 
which (taken over all distinct symbol vector pairs {s, s'}) 
must be maximized. 

At high SNR, the upper bound in Eq. iT% is dominated by the 
exponent —rMn of p. This leads to a more general definition 
of diversity gain [33] as being the total degrees of freedom 
offered by a communication system, reflected by the factor 
involving the negative power of the SNR in the expression 
of the error probability. Full diversity gain is achieved when 
the total degrees of freedom {— MM ft) offered in the multi- 
antenna system are utilized. We adopt this latter notion of 
diversity gain when we consider the STBC for the MISO 
system. 

Since ilZ/j = 1, full diversity for a MISO system is achieved 
if the exponent of the SNR in the expression of the error 
probability is equal to —M. Let us now consider the condition 
on Ti, for which full-diversity is achieved by a MISO system 
employing a linear receiver We need only to consider the 
use of a linear zero-forcing (ZF) receiver because the same 
condition extends to MISO systems using linear minimum 
mean square (MMSE) receivers or other more sophisticated 
receivers. Since the diversity gain of a communication system 
relates the probability of error to SNR, we first analyze the 
symbol error probability (SEP) of detecting different signal 
constellations by a linear ZF equalizer and express these in 
terms of the SNR. 

A. Symbol Error Probability of Various Signalling Schemes 

Here, we examine three commonly used signalling schemes: 
1) square QAM, 2) PAM and 3) PSK constellations respec- 
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lively. Let /i denote the cardinality. Firstly, we summarize 
the definition of some common parameters which govern the 
performance of the ZF linear detectors under these schemes. 
We use the index i = 1, 2, 3 to denote parameters associated 
with the three signalling schemes as ordered above. Let 
Esi, i = 1, 2, 3, denote the respective average symbol energy 
in each of the above schemes, and let a^ be the noise variance 
at the receiver antenna. Therefore, the SNR for each symbol 
at the receiver is given by 

p, = E.,/a^; i^l, 2, 3 (13) 

Note that a [Ti. 'W^ is the noise power at the output of the 
ZF equalizer for the £th symbol. 

1. Square QAM signals: The SEP of a ZF receiver for the 
square QAM signal s^ is [46] 

Pi(h,S£) = 



2. PAM signals: The SEP of the ZF receiver for a /^(-ary PAM 
signal si is given by [46] 



P2(h,s,) 



Now, using Eq. ( fTSl l and noting that by putting sin 6* — 1 
in the integral, we have Q{z) < exp(— z^) for z > 0. 
Together with p2 = Es2/<^'^, we arrive at an upper bound 
of Eq. (O, 



/ 








\ 


2(M-1) 


.\ 


iE,2 


(^2- 


-l)cr2 


( 


u 1 
19) 



i^2(h,Sf) < ^ exp 




3p2 

(20) 
3. PSK signals: The SEP of a ZF receiver for the PSK signal 
SM is given by [46] 



P3(h,s^) = - / exp 



( 



EszS\t?{'k/p) 



\ 



2a2 



n^H 



sin^ e 



de 



where Q{z) — ^== J^°° e ^ ^'^dx. We use the following 
alternative expressions for the Q- and Q^. functions [46]- 
[49] 



/ 
(21) 

Writing p2 = Es2/<y'^, similar to the PAM signal, Eq. (1211 ) 
can be upper bounded by 





1 


.71/2 


Qiz) = 




/ exp 




■K 


h 




1 


r^/i 


r(z) = 




/ exp 




■K 


h 



2 sin^ e 



P3(h,s,)<^^^^exp 



d9 z > (15) 
de z > (16) 



P3sin^(7r//x) 



\ 



n"H 



(22) 



J 



2sin2 6' 

Substituting Eqs. (fTSb and ( fT6] ) into Eq. ( fT4l i and after a 
little manipulation, we obtain 

4 / 1 



Pi(h,Sf) 



1 - 



^/M 



7r/4 

exp 



/ 







3i?si 



\ 



/ 

Tr/2 

exp 

7r/4 



1 - 



4(p-l)fT2 
1 



(-H^T^) ' 



sin 



u 



iEsi 



\ 



V 



4(^-1)0-2 



{h"u) ' 



sm 



^6* (17) 



B. Design Criterion for Full-Diversity STBC for a MISO 
System with Linear Receivers 

We now examine the diversity gain achievable by a MISO 
system. 

Theorem 1: For a MISO system employing a square QAM, 
a PAM, or a PSK signalling scheme of cardinality fi in 
the transmission, a linear receiver (ZF/MMSE) achieves full 
^Q diversity for the system if 7i Ti. is non-singular for any 
nonzero h, or equivalently, if X^(s)X(s) is non-singular for 
any nonzero s. ■ 

Proof: From Eqs. ( fTSl l, ( l20t . and ( |22] ). we can arrive at a 
generalized upper bound on the symbol error probability for 
the /i-ary QAM, PAM, and PSK signals such that 



ee 



7 



Pi{ii,si) < exp 



at Pi 



We can obtain an upper bound for Eq. ( [TtI i by putting 
sin6' = 1 in the two integrals and, writing pi = Esxja^, 
this easily simplifies to 



71^71 



J 



1, 2, 3 (23) 



Pi(h,s,)< 



M 



M 



■ exp 



3pi 



4(m-1) 



Ti^Ti 



(18) 



where p^ — Esi/cr^, and 

fli = 3/[4(/x-l)],a2 - 3/[2(/^2_^)]^^j^^^^ ^ sin2(Vp)/2 

(24) 



jH 



Since ?i ?i is non-singular for any nonzero h, we can apply 
Eq. (|9]) on Eq. ( |23] ). Here, we see that the arithmetic mean 
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of the SEP of all the three signalling schemes have a general 
upper bound given by 



Pi{h) < ^^^exp(-a,p,Co||h||2) 



M 

A^-1 
M 



cxp (-aiPiCoh^h) 



(25) 



Now, h is assumed to be Gaussian with zero mean and 
covariance matrix S. Therefore, averaging the exponential part 
of the right side of Eq. dZSl l over the density function of h 
yields 

^MdetS / """P {-^^P^Coh"h) cxp {-h^'^-'h) dh 

^^ "^ ^^ = dot I + a,;p,CoSr ' (26) 

Substituting Eq. ( l26b into Eq. dZSl l. we establish the following 
inequalities: 

M-1 



E[P.(h)] 



< 



< 



A^-1 



■dct(I + ajp,CoS)-i 



det(CoS)-iar''' ) P,-''' 



z = l,2,3 (27) 

The exponent of pi in Eq. ( |27] l indicates that the upper bound 
of the SEP using a ZF receiver in a MISO system to detect 
signals from the three schemes indeed achieves full diversity 
for non-singular 7i 7i. 
We now show the equivalency of the following two statements; 

1) Ti Ti is non-singular for any nonzero h; and 

2) X^(s)X(s) is non-singular for any nonzero s. 

We will show 1) =^ 2), and the reverse can be similarly proved. 
From the development of Eq. (|5]l, we have X(s)h = Tis. 
Now, if 7i 7t is non-singular for any nonzero h, then 7i 
has full column rank, and hence Tis 7^ V s 7^ 0. Therefore 
X(s)h 7^ for any s 7^ 0, h 7^ 0. This implies full column 
rank of matrix X(s), and hence X^(s)X(s) is non-singular 
for any nonzero s. D 

Remarks on Theorem 1 : 

a) Although the proof provided here is for square QAM, PAM 
and PSK signallings. Theorem [T] can be shown to be valid 
for any signal constellation. 

b) Since the condition provided here is sufficient for a linear 
receiver to achieve full diversity, the same condition natu- 
rally yields full diversity for more sophisticated receivers 
such as MMSE/ZF-DFE or ML receivers. 

c) For a MIMO system using an ML detector, the requirement 
for full diversity as indicated by Eq. (fT2l) is that the coding 
matrix X^(e)X(e) is maintained at full-rank for the 
signals s, s' e S. However, for a MISO system employing 
a linear receiver. Theorem [T] shows us that full diversity is 
achieved if the coding matrix X^ (s)X(s) is of full rank for 
any signal s, a much stronger condition than that required 
by systems using an ML detector 

In the following section, we present the Toeplitz STBC which 
has a simple structure satisfying the full-rank condition in 
Theorem [U and is therefore a full-diversity STBC for a MISO 
system employing a linear receiver. 



IV. Toeplitz Space-Time Block Codes and Their 
Properties 

A. Toeplitz STBC for a MISO System [41] 

To examine the structure of the Toeplitz space-time block 
code, we let a. — [ai a2 ■ ■ ■ cxl]'^ ■ A {K + L — 1) x 
K Toeplitz matrix generated by a and a positive integer K, 
denoted by T{a., L, K), is defined as 



[T{(x,L,K)],,^ 



0, 



if j > j and i — j < L 



which can be explicitly written as 

/ ai ... 
a2 Oil 



T{a,L,K) 



a.2 







otherwise 



\ 




"1 

a2 



(28) 



V 



(29) 



aL / 



{K+L-l)xK 



If we replace a by s, the information symbols to be transmit- 
ted, then a Toeplitz STBC matrix A'b(s) is defined as 



A'B(s) = r(s,L,X)-B 



(30) 



where, for K < M, B is a K x M matrix of rank K placed 
in the coding matrix to facilitate the transmitter antennas with 
beamforming capability. At time slot n, the ?ith row of the 
N X M matrix Ab(s) is fed to the M transmitter antennas 
for transmission. Apply the Toeplitz space-time block coding 
matrix to the MISO system described in Eq. (O, and we have 



y = T(h,/i,L)s + ^ 



(31) 



where h = Bh, and L = iV - iiT + 1. Thus, T{h, K, L) can 
be viewed as the overall channel matrix of the MISO system. 
Example 2. For K = M = L = 2,N = K + L - 1 = 3, 
and B = I2, the codeword matrix and channel matrix are, 
respectively. 



^-1,(8) 




r(h,2,2) 




For this code, there are L = 2 symbols to be transmitted in 
A^ = 3 channel uses. Therefore, the symbol transmission rate 
of this system is i?^ = | symbols per channel use. 
Remarks on Toeplitz STBC: 

a) Eq. ( I3TI ) is identical in form to that describing a MIMO 
intersymbol interference channel for zero-padding block 
data transmission (e.g. [50]). It can thus be interpreted 
that the original MISO channel is transformed into a 
Toeplitz virtual MIMO channel. In other words, the space 
diversity has been exchanged for delay (time) diversity. 
This is reaUzed by transforming the flat fading channel 
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into a frequency selective channel with zero-padding. This 
technique is parallel to that employed in [40].!^ 

b) For such a system, we can utilize the efficient Viterbi 
algorithm [51] to detect the signal s if perfect channel 
knowledge is available at the receiver. On the other hand, 
when channel coefficients are not known at the receiver, we 
can make use of the second order statistics of the received 
signal to blindly identify the channel [50], [52]. 

c) Toeplitz STBC is a non-orthogonal STBC whose coding 
matrix X^X possesses non-vanishing determinant for any 
signalling scheme. Hence according to Theorem[T] the code 
achieves full diversity even with the use of a linear receiver. 
On the other hand, since full diversity STBC designed 
for ML receivers (e.g., [34]-[36], [53]) maintain non- 
vanishing determinant only for certain types of signalling, 
full diversity gain is not guaranteed when a linear receiver 
is used. 

d) When B = I, Toeplitz STBC becomes a special delay 
diversity code (DDC) [37], [38] with padded zeroes. In 
general, DDC is applied with the use of outer channel 
coding and ML detectors to achieve the full diversity gain. 
However, here we show that the Toeplitz STBC possesses 
special properties which enable full space diversity to be 
achieved even with the use of the simplest Unear receiver 
and the signals can be of any type. 

B. Properties of Toeplitz STBC 

We now examine some important properties of the Toeplitz 
space-time block codes introduced in the previous subsection. 
These properties will be useful in performance analysis and 
code designs in the ensuing sections. 

Property 3: The definition of the Toeplitz space-time code 
shows that the symbol transmission rate is i?^ = -^ = '^^l^^^ 
symbols per channel use when K < M. Therefore, for a fixed 
AI, the transmission rate R can approach unity if the number 
of channel uses is sufficiently large. 

Property 4: For any nonzero vector a, there exists 



< Ct, 



< 



Cti 



< 1, and the matrix 



(T^ {a, L, K)T{a., L, K)) satisfies the following inequality. 



Crminllaf^ < dot {T" {a, L,K)T{a,L,K)) 



< Ct. 



\2K 



(32) 



Proof: By letting a = h in Eq. Q and choosing 



0(L-l)xl 
I(L-l) 



where 

Lllx(L-l)/ \U(A'-l)x 

we obtain an equivalent channel T-t of the same structure as 
T{a., L, K). Hence, T{a., L, K) is a special case of T-L. Thus, 
from Property [T] there exist Ct min and Ct max for which 
Eq. (|32] | holds. Now, we note that the diagonal entries of the 
matrix T^ [a, L, K)T{a., L, K) are all the same and are equal 

'in [40], consideration is given only to the use of an ML detector and 
performance analysis is restricted to BPSK signalling under the assumption 
that only one bit error occurs. 



to [T"{a,L,K)T{a,L,K) 



kk 



a 



fc = 1, • • • , A' . By 



applying Hadamard's inequality [44], we arrive at: 

dot {T"{a.,L,K)T{cy.,L,K)) < Ijaf^ 



Cti 



< 1 



(33) 



Furthermore, since a is nonzero, we can assume, without loss 
of generality, that the first element ai 7^ 0. (Otherwise, we 
can always permute the nonzero element to the first position.) 
The N X K "tall" matrix, T{a, L, K), can be partitioned into 
a top K X K matrix J7i, and a bottom matrix ^2 containing 
the rest of T{a., L, K), i.e., 

/ ai ... \ 



T(a, L, K) = 



a-i 



an 



ai 
an -I 



"A'+l 



OiL 





aL-1 







0-2 

CHL-K-l 
OiL-K 



Oi 



O2 



V ••• aL J 

We note that fii is a lower triangular matrix having equal 
diagonal elements qi 7^ 0. (Here, we assume that K < L. The 
proof is equally valid if i < K by exchanging the roles of K 
and L). Since T" {a, L, K)T{a, L, K) = Oi^r2i +02-^^2, 
using a standard result [44, p. 484] on the determinant of the 
sum of a positive definite and a positive semi-definite matrix 
we have 

det (T^(a, L, K)r{a., L, K)) 

> det(f2i^r2i) + dct(r22^f^2) 

> dct(f2i^r2i) 

= |«i|''^ (34) 

Since cvi 7^ 0, dct (T^{a, L, K)T{a, L, A')) is nonzero for 
any nonzero a and thus we have 

Crrain > (35) 

The conclusions in Eqs. ( |33] | and (|35] | complete the proof of 
the property. D 

We see that T is non-singular by Property |4] then using the 
fact that T as a special case of H in Property |2] we have: 
Property 5: 

[{T"{a,L,K)T{a,L,K)yXk ^ C^T„nn||a||2 

for fc== 1,2,--- ,A' (36) 

where the existence of Ct min > is shown in Property |4] 
holds. ■ 

We now introduce the following definition related to the 
measure in a signal constellation S: 

Definition 1: For s,s' G S, the minimum distance of the 
signal constellation is defined as 



,{S) = mill I 



(37) 
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If II s — s'll = d„iin{S), we say that s and s' are neighbours. 

From the definition of the coding matrix X in Eq. dSOl l with 
B = Im, we can now establish a lower bound for a metric 
between Aij,j(s) and A'ij^(s') in dnuniS). Let e = (s - s'). 
For notational convenience, we let Xij^j{e,{ii,i2,- ■ ■ ,im}) 
denote the matrix consisting of m of the columns of Aj^,^ (e) 
indexed by {ii, Z2, • • • , im} where ii < i2 <,■■■,< i,n and 
these columns are not necessarily consecutively chosen. Then, 
we have 

Property 6: For s 7^ s' e S^, where S^ = S x S ■ ■ ■ x S, 
and any nonzero vector e = (s — s'), we have 



det [A'i^^(e, {ii,i2,- ■ ■ ,ir,i})Xi^,{e,{ii,i2,- 



.})] 



><-(5) (38) 

for m = 0,l,-- , M —1, where the equality holds if and only 
if s and s' are neighbours, i.e., iff ||e|| = (i,„in(5). ■ 

Proof: The proof of this property is similar to that of Prop- 
erty |4] As in Property ID without loss of generality, we can 
always assume that ei 7^ with ei being the first element of 
e. Thus, Aij^j (e) can be written as 



62 



SM 



XuA^) 





ei 

62 



gm 







\ 



V 







ei 

62 

en 

6L 



(39) 



/ 



NxM 



An important observation in Eq. ( I39b is that the top submatrix 
consisting of the first M rows of Xij^^ (e) is a 7\f x M 
lower triangular matrix with nonzero diagonal entries and 
therefore, nonzero determinant. We can also see that the sub- 
matrix Xij^, (e, {Ji, J2, ■ • ■ , im}) preserves the same property 
because by permuting its rows and columns, an rn x m lower 
triangular matrix can always be formed as its top part, i.e., 
A'ijj(e, {ii, Z2, • • • , im}) can be expressed as 



niA'i„(e, {ii,Z2, 



5 ^r 



J)n 



2 — 



I e^ 



where Hi and 112 denote the N x N and m x 
m permutation matrices, respectively, Jli contains the 
first m rows of Aij,^(e, {11,^2, •• • ,im}) and hence is 
lower triangular, and 02 denotes the remaining subma- 
trix of Xi^{e,{ii,i2,- ■ ■ ,im})- Since the permutation of 
the rows and columns of Xij^j{e,{ii,i2,- ■ ■ ,im}) does 
not change the determinant of its autocorrelation ma- 
trix Xf^^{e,{ii,i2,--- ,im})Xi^i{e,{ii,i2,--- ,im}), there- 
fore, as in Property IH we arrive at 

'ist ['^iM(e,{ii,«2,--- ,im})Xi^,{e,{ii,i2,--- ,im})] 
> det(f2i^f2i) + det(f^2''^2) > <1„(5) 



where the equality holds iff ^1,2 is a zero matrix, i.e., iff s and 
s' are neighbouring points. D 

We can establish another useful property on the metric 
between Aj^^^ (s) and Aj^,^ (s') by first recalling an important 
property in matrix algebra [54]: 

The characteristic polynomial of an M x M matrix A is the 
polynomial whose roots are the eigenvalues of A. Mathemat- 
ically, it can be re- written as 



h{v) = dct (I + i/A) 



sucli that 



v^' + c,v^'-' 



CM-lV + CM 

^det(A),,,. 



(40) 
(41) 



where A, 



denotes the principal submatrix obtained by 



deleting the rows and columns of A except the zith, the Z2th, 
• • • , and the imth ones, and -d denotes the combination set 
of ii, • • • , im- We note, in particular, ci = tr (A) and cm ~ 
det (A). Now, the following property provides us with another 
lower bound on the metric between Ai„(s) and Ai„(s') in 
relation to dmin(5). 

Property 7: Let A = diag((5i, (^2, • • • ,(5j\f) with 5„i > 
for m = 1,2, •• • ,M. Then, for any nonzero vector e, the 
following inequality holds 

M 

det (A + X,l {e)Xj,, (e)) > J] (-5™ + <in(5)) (42) 

m— 1 

with equality holding if and only if s and s' are neighbours. ■ 
Proof: Let us first rewrite the left side of Eq. ( |42| ) as 



det(A + <^(e)Ai,,(e)) 



det(A)det 1+ A'i,,(e)A-i/2 Ai,,(e)A-i/2 (43) 



H 



.H 



Now, let A = (A'i^,(e)A-i/2)" A'i^^(e)A-i/2, then 
Eqs. ( |40l i and W\\ becomes 



det 



I 



H 



Ai„(e)A-i/2 Au,(e)A-i 



/2 



= 1 



M 
?n — 1 



where 



Cm = ^det [Ai^,(e, {ii, 



«2,--- 



J 1-m,s) 



A'i„(e,{ii,i2,--- ,«m})]n^»/ (45) 
(=1 

Using Eq. dJSl ) in Property |6] on the right side of Eq. ( |45] |. we 
have 



>d'^^{S)Y,l[5-^ 



(46) 



1? e=i 



where 'J2§ Y[T=i ^i~^ denotes the sum of the combination of 
the product of (5^^ taken to at a time. Equality in Eq. ( |46] | holds 
if and only if s and s' are neighbours. Combining Eqs. ( l44b 
and ( |46] | results Jn 



det 



1+ A'i,,(e)A-i/2 A'i,,(e)A-i 



H 



>-l/2 
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> 1 



M 



^d^::in{s)Y.X{K' 



■0 £=1 



m— 1 









(5)) 



(47a) 



(47b) 



where in the second step, we recognize that the right side 
of Eq. ( 147 at is the eigen polynomial of the matrix d^^^jj^A 
which, in turn, equals det(I + d^^^A.) and hence Eq. ( |47b| ). 
Combining Eqs. ( |47b| ) and ( |43] l, we complete the proof of 
Property H D 

V. ToEPLiTZ STBC Applied to a MISO System with a 
Linear Receiver 

We now apply the Toeplitz STBC to the MISO communi- 
cation system using the properties presented in Section IIV-BI 
From Property [3] and |4] in Section IIV-BI and Theorem [T] in 
Section ITlI-B I we can see that the Toeplitz STBC can approach 
unit-rate as well as full diversity even if only a linear ZF 
or linear MMSE receiver is used in a MISO system. In the 
following, we examine the optimal tradeoff between diversity 
gain and multiplexing gain [33] when the Toeplitz STBC is 
employed in a MISO system equipped with a linear receiver 
We first make the assumption that the channel coefficients are 
independent, i.e., S = I. 

Now, our MISO system has M transmitter antennas trans- 
mitting a signal vector s of length LinN = L + AI — 1 
time slots. Also, all the above three signalling schemes have 
constellation cardinality fi. Thus, employing any of the three 
schemes described in SectionlUin our MISO system will result 
in a transmission data rate t given by 

r = ^ logs M (48) 

Note that r is the bit rate of transmission. The multiplexing 

gain g, on the other hand, is dependent on the scheme and in 

general, is defined as [33] ^ 

.9. = , ^^T^ (49) 

log2 bNRj 

where "SNRi" refers to the general SNR in the received data. 
Here in our analysis, we use the SNR of the received data 
block for the SNR^ and denote this by pbii, ^ = 1, 2, 3 when 
the ith signalling scheme is employed. Notice that in the MISO 
system, we always have < g,; < 1 V i, since the system has 
only one receiver antenna. Hence, we can write /i = p^if ■ 
Thus, if we want to maintain nonzero multiplexing gain at high 
SNR, then fi, the cardinality of the constellation must be large. 
It has been shown [33] that at high SNR, we can trade-off the 
multiplexing gain for diversity gain and vice versa, and the 
optimum trade-off for our MISO system with M transmitter 
antennas is given by 



Dop = M{1 - g) 



(50) 



where Uop is the optimal diversity gain. Let us examine trade- 
off of the multiplexing gain for diversity gain in the three 
signalling schemes: 

1. QAM signals: The averaged symbol energy Eg for square 
QAM signal is [55] 

Esi^\Ui-l) (51) 



We note that Egi increases with the constellation cardinal- 
ity ^. From Eq. ( fSTI i, the averaged transmission energy per 
block can be calculated as E'bu = |(/i — \)ML. Given 
a^ being the noise variance at the receiver antenna, the 
averaged noise power per block is cr^jj = a'^N . Therefore, 
the block SNR is 

2(/i - l)ML 



Ph\i 



iNa^ 



leading to 



2 2nML 2ML wai 1 
O ~ -rrz = -TT^Pbll ^ 



(52) 



STVpbu iN 
where the approximation is under the assumption of large 
/i, and Eq. ( |49] l has been used. Therefore, from Eqs. ( fT3l l, 
(l24b , dSTT ) and (|52]i, we obtain 

?,Es 1 3iV 



flipi 



PhM 



(53) 



4(/i - l)cr2 20-2 AML 

Now, consider Eq. dZTJ i on the upper bound of the SEP for 
a ZF receiver, i.e., 

E[Pi(h)]<Cj;^^.„(aiPi)-*^ (54) 



(55) 



Substituting Eq. (I53t into Eq. (I54b . we obtain 



-M 



E[P.(h)l<c-4,|^^) «.-'-' 



Hence, the diversity gain for this scheme, Di, is given by 



Di{g)^M[l 



N \ 
- -gj = M{1 - g) - eMg 

= Dop{g)-eMg (56) 

> 0. From Eq. (|56T l, we can see that 
Dop{g). However, we can make e small by 



M-l 



where e = 

^1(5) < 

choosing L sufficiently large so that Di{g) w Dop{g)- 
Hence, Di{g) is the e-approximation of Dop{g)- We can 
always choose L = [ ^^~"'" ] + 1, where [•] denotes the 
integer part of a quantity, and therefore, we can say that 
the ZF receiver is able to approach the optimal diversity- 
multiplexing tradeoff if the proposed Toeplitz code is used 
with a square QAM signalling scheme of large cardinality. 
If a MMSE receiver is employed, utilizing an expression 
parallel to Eq. dSST l for the MMSE receiver, we can also 
show that the optimal diversity-multiplexing tradeoff can 
be asymptotically achieved if the Toeplitz STBC is applied 
to the MISO system in which a square QAM signalling 
scheme of large cardinality is used for transmission. 
PAM signals: We note that the averaged transmission 
energy Eg for /i-ary PAM signal is given by [55] Eg = 
i(/x2 — 1). Hence the averaged transmission energy per 
block is Eshi = |(m^ ^ 1)ML. Following similar argu- 
ments resulting in Eq. (|52] |. for PAM signals we have. 



a « — 

6A^Pbi2 

Also, similarly to Eq 
ZEs 

0.292 



2(^2 _ 1)^2 



hlL 2]V3 

= -^Pbl2 - 

, we have 

1 _ 3iV 
4^^ 



2ML 



Ph\2 
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Therefore, from Eq. jTTI i. the upper bound on SEP for PAM 
signal is 



E[P2(h)] < C^:L(a2P2) 



-M 



^ ( 



3NCt 



\ 2ML 
Hence, the diversity order is 

2N 



-M 



Pbl2 



-g-M 



D2{g) = M ( 1 - —g\ = M{1 - 2g) - 2eMg 

< A/(l - 2g) < M(l -g) = D,p{g) (57) 

EquaUty in dSTJ i holds iff g = 0. Therefore, for finite 
multiplexing gain, D2{g) cannot approach the optimal 
tradeoff Dopig). 
3. PSK signals: The averaged transmission energy Eg for 
fi-nry PSK signal is [55] Eg = 1. Hence the averaged 
transmission energy per block is Esu = ML. Therefore, 
we have 

ML 



Nphl3 



We also have 



0-3 P3 



Eg sin {tt/ p.) 
2^ 



Nn' 



1 2Ng 

PbVi'-— (58) 



2cr2^2 2ML 

where, the second step comes from the assumption of large 
/x. Following similar arguments as PAM scheme, it can 
be shown that PSK signalling cannot achieve the optimal 
tradeoff of diversity-multiplexing gains. 

VI. Optimal Toeplitz STBC Design for MISO 
System with ML Detector 

The previous section shows what could be achieved when 
the Toeplitz STBC is applied to a MISO system equipped 
with a linear receiver. In this section, we will examine the 
application of the Toeplitz STBC to a MISO system equipped 
with a ML detector. In particular, we seek for the optimal 
design of the matrix B inherent in a Toeplitz space-time 
block code Eq. (|30] | such that the worst case pair-wise error 
probability is minimized when a maximum likelihood detector 
is employed. 

Given a channel realization h and a transmission matrix B, 
the probability P (s -^ s'|h, B) of transmitting s and deciding 
in favor of s' 7^ s with the ML detector is given by [56] 



P(s^s'|h,B) = 



ff(s,s') 
2ct 



(59) 



where d{s,s') is the Euclidean distance between the Toeplitz 
coded signals s and s' after being transmitted through the 
channel, i.e., it is the Euclidean distance between ^3(8)11 and 
A'B(s')h. Because of the relation of Eq. dSOl l. we can write: 



d^(s,s') = (s-s')''r^(h,M,7V)r(h,M,7V)(s-s') 
= h"Xi{e)XB{e)h (60) 

where e = s — s'. By employing the alternative expression of 
the Q-function in Eq. (flST l and taking the average of Eq. (|59] l 



over the Gaussian random vector h, the average pair-wise error 
probability can be written as 



P(s^s'|B) = 



1 



7r/2 



d9 



det I ■ 



ia"^ sin^ ( 



^SA'^(e)A'B(e)) 



(61) 



with S being the covariance matrix of h. Our design problem 
can now be stated as: 

Design Problem: For a fixed number M of transmitter an- 
tennas, find a K x M, {K < M), matrix B such that the 
worst-case average pair-wise error probability P(s -^ s'|B) 
is minimized, subject to the transmission power constraint, 
tr (B^B) < 1, i.e.. 



B 



op 



arg mm 

tr(B"B)<l 



max 

1, s'es' 



P(s^s'|B) 



(62) 



where S'" = S x S ■ ■ ■ x S with S denoting the signal 

constellation of each element of s. 

To solve the above design problem, we not only have to 

find the optimum B, but also have to determine its dimension 

K. Let us first examine the M x M covariance matrix, S = 

E[hh ], of the transmission channels h. Suppose we perform 

an eigen decomposition such that S = VA V^ where V is an 

M X M unitary matrix and A = diag(Ai, A2, • • • , Xm) with 

Ai > A2 > • • • > Aj\/ > 0. The following theorem provides 

us vi-ith an Qptirnuni_,desigjj of B: , ^^ ,^ , 

Theorem 2: Let T = t[Tag(7i, 72, • • • ,7_ft:j, A < Ai, be 

the singular values of B and let G{A.kT, e) denote the integral 

GfA^r, e) = - / TT I 1 + ^-^ I ^6* for e > (63) 
^ J^ fc=i V sin^ 9 J 



We obtain an optimal T by solving the following convex 
optimization problem^ 



:'op = arg min G A^F, ■ 

tr(r)<l ' 



.{S) 



8a' 



(64) 



where Fop is a A' x K diagonal matrix given by Fop = 
diag(7opi, 7op2, ••• , 7opA') with K being the highest 
integer for which [Fop]fcfc = 7opfe > 0, fc = 1, 2, • • • , AT. Then 
the optimum transmission matrix is given by 

(65) 



Bop — A op » x 



where V/^ is the MxK matrix containing the K eigenvectors 
corresponding to the K largest eigenvalues in the eigen 
decomposition of S. Furthermore, the worst case pair- wise 
error probability is lower bounded by 



max P(s 



s'|B) > G AifFopt 



'^min('^) 
8^2 



(66) 



- Note that the work presented here is different from that in [57] in which 
a precoder matrix is designed for a frequency-selective fading channel even 
though both involve Toeplitz structured matrices. Here, the Toeplitz matrix 
containing B and h is separated from the signal vector s. This, by the 
properties of the Toeplitz STBC shown, transforms the design of B into 
a convex optimization problem. In [57], however, the design parameter and 
the signal vector are all parts of the Toeplitz structure resulting in a non- 
convex design problem that can only be solved by numerical method with no 
guarantee for global optimality. 
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Equality in Eq. ( |66] l holds if and only if 

i) ||s-s'|| = dminiS), and 

ii) B = Bop. ■ 

Proof: We first establish a lower bound on the worst case 
average pair-wise error probability. Let s^ and s^ be neighbour 
symbol vectors differing in only the ^th symbol, i.e., S£ — s^ = 
e£ = [0 • • • eg ••• 0]-^ where \eg\ = dniin{S). Then, 
we can write 

8cr^ sm ff 



M 



< 



det(lM 

fc=i 



mill V 



80-2 sin-" 9 

'^min('5) 



-SB^B 



8cr2 sin^ e 



Ykh 



(67) 



where the first step is a result of the structure of eg on 
the Toeplitz code, and the second step is the result of an 
inequality for the determinant of a matrix [58], [59]. Equality 
in Eq. (E) holds if and only if B = B^ = TVa-^, i.e., the 
singular vectors of B are the eigenvectors of S. Substituting 
the inequality of ( |67] i in Eq. (|6Tl i. we have P{si -^ s^'jE) > 



Since 



max P(s -^ s'|B) ) > 



P{s( -^ s^'|B), the worst case average pair-wise error 
probability is lower bounded by 



max P(s ^ s'|B) > G (akT, "^"f f ^ 



(68) 



If we minimize both sides of Eq. ( I68l l, we can write 
minf max P(s ^ s'lB) ) > G | AxFop '^"""^'^^ 

B \s,s'e5^,s^s' / \ 



8a' 



(69) 
where Top is obtained according to Eq. (|64l i. 

Let us now establish an upper bound for the worst case 
average pair-wise error probability for the specially structured 
transmission matrix Bq above. For any error vector e, we have 



det Im+ 



1 



8cr2 sij^^ 



SA'#fe)A'B„(e) 



1 



80-2 sin-" e 



M 



dciAKT')dct{A + X,l{e)X,,,{e)) 



(70) 



where the special structure of Bq has been utilized, Ak 
denotes the diagonal matrix containing the largest K positive 
eigenvalues of S and A = (8cr2 sin^ 6')A^^r"2. Using 
Eq. ( iTOt in Property |7] for any nonzero vector e and nonzero 
6 in the interval [0, 7r/2], we have 



M 



det(A + <^(e)A-i„(e))>n 



fe=i 



BfT^ sin^ ( 



7^A. 



.(5)) (71) 



where, according to Property |7] equality holds if and only if 
s and s' are neighbour vectors. Eq. ( iTOt and Eq. ( fTTT i together 
yield 



detd 



1 



M 



8cr2 sin2 



SA'#Je)A'B„(e)) 



> 



k=l 



8cr2sin2 6' ' 



(72) 



Again, substituting Eq. ( |72] ) in Eq. jM] ) and using the optimum 
Top yields 

max P(s ^ s'IBop) < G (akT^p, %i^^ 

where equality holds if and only if ||s — s'j| = (imin('5). This 
results in 



mill I max P(s -^ s'|B) 
< max P(s -^ s'jBop) 

<r r ( \ r ^min(^) 



Combining Eq. ( l73T l with Eq. ( |69] l yields 



s'lB) =G A^r 



k-XJ-op, 



.(5) 



(73) 



(74) 



mm I max P(s 

Eq. dllli holds iff B = TopVf and ||s-s'|| = d,„i„(5). Thus, 
the proof of Theorem |2] is complete. D 

Remarks on Theorem El 

a) Theorem|2] shows us that the lower bound of the worst case 
pair-wise error probability can be reached by having B = 
FopVl^. Thus, the design problem in Eq. ( |62l ) becomes 
finding the optimum Fop of Eq. ( l64l i. 

b) The original non-convex optimization problem has been 
transformed into the convex problem in Eq. ( |64b and can be 
solved efficiently by interior point methods. The convexity 
of the objective function can be verified by re-writing 
G(AKr, e) in Eq. ^ as 



G{AKT,e) 



We notice that — In(-) is a convex function over 72, and 
hence their sum is also convex over [71 , • • • ,7^]. Now, 
exp(a;) is monotonically increasing with x. By composition 
rule [60] (Page 84), the integrand in Eq. jTST l is a convex 
function implying that G(A/^r, e) is convex, 
c) The solution of Eq. ( l64l i yields the values of the diagonal 
elements {7opi, 7op2, • • • , 7opA/}. Some of these values 
may not be positive. We choose all the K positive ones to 
form the singular values of Bop. 
Theorem 12] provides us with an efficient scheme to obtain the 
optimal matrix Bop by numerically minimizing G(A/i-r,e). 
However, if the Chernoff bound [46] of the pairwise error 
probability is employed as the objective function for mini- 
mization instead, a closed-form optimal B can be obtained. 
This can be shown by setting sin^ = 1 in the pairwise error 
probability of Eq. dMT l, so that we obtain the Chernoff bound 
as 

P (S ^ S'|B) < „,,, ,, , ,„ ,^_L..„._^.. ._.N (76) 




2dct (I + (8cr2)-iSA'^(e)A'B(e)) 
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Seeking to minimize the worst case Chernoff bound, and fol- 
lowing similar arguments which establish the the optimization 
problem in Eq. ( |64l ). we arrive at the following problem, 



M=4, L=8, Different Signal Consteilation 



■ op 



arg rnm 

tr(r)<i 2 



1 ^' 



fe=i 



AS) 



>^kil 



(77) 



where Fop is a diagonal matrix with diagonal elements 7opfe. 
This problem is a relaxed form of that in Eq. (|64] | and its 
solution is provided by the following corollary: 

Corollary 1: The solution. Fop, for the optimization prob- 
lem of Eq. ( ITTI i can be obtained by employing the water-filling 
strategy [59]. The diagonal elements of Fop are given by 



Topfc 



8a2 



\ M 



Ma 



dL^iS) ^ Xi 



Mo , 

El, 



1 

1, 



-I + 



,K (78) 



where notation [a;]+ denotes max(a;,0). The optimal choice 
of K is K ^ Mq, where Mq is the maximum positive integer 
satisfying 




8a^ ^l\ 1 



m = l,2,--- ,Mo 



The optimum transmission matrix is thus Bop — FopV^. 
Proof: The minimization of the type of problems of Eq. dTTj i 
has been studied by several researchers and the water-filling 
solution is shown in [59]. D 

Remarks on Corollary [T] 

a) For the particular case in which the channel coefficients 
are mutually independent, i.e., S = Im, then any AI x M 
unitary matrix scaled by a factor l/\//\7 is a suitable choice 
for Bop. 

b) The optimal design Bop maximizes the coding gain [45] 
which is defined to be the normalized minimum determi- 
nant of X{e)^ X{e) for all nonzero e. The criterion of 
coding gain maximization is derived from minimizing the 
Chernoff bound on pairwise error probability for a system 
equipped with a ML detector. Since Bop minimizes the 
worst case Chernoff bound on PEP, we conclude that it 
achieves optimal coding gain. 

Remarks on the Optimum Transmission Matrix Design: 
• The derivations of Theorem |2] and Corollary [T] are based 
on the consideration of using a ML receiver in the MISO 
system. For the case when the system is equipped with a 
linear ZF (or MMSE) receiver under the environment of 
correlated channels, the problem of obtaining an optimum 
B becomes very complicated. This is because we seek 
for a matrix B to minimize the respective average error 
probability obtained by averaging the expressions of error 
probability in Eqs. ( fT4b . ( fT9l ) and ( l2Tl i for the respective 
signalling schemes over the random channel matrix. This 
requires the knowledge of the PDF of the equivalent 
channel matrix. However, the equivalent channel matrix 
in these cases is of a Toeplitz structure for which the 
PDF is unavailable, and therefore, the expressions for 



4-QAM 
16-OAM 
64-QAM 




SNR [dB]/bit 



Fig. 1 . The average BER performance of the proposed ToepUtz STBC when 
signals are selected from different constellations. 

the average error probability cannot be obtained. (For the 
case of linear MMSE receivers, a similar problem exists). 
• An alternative way to attack the problem in the case when 
a linear receiver is employed is to consider the upper 
bound on the averaged error probability given in Eq. ( l27b . 
We can minimize this upper bound with respect to B. 
However, this necessitates the knowledge of the value 
of Co- From Property H] Co is the minimum value of 
the determinant of the Toeplitz matrix having its column 
vector belonging to unit ball and this renders the solving 
of Co difficult. Thus, in this paper, we are unable to come 
up with any true optimal B for MISO systems equipped 
with linear receivers. 

VII. Numerical Experiments 

In this section, we examine the performance of the Toeplitz 
code in a MISO system. We first evaluate the performance of 
the system equipped with a linear receiver under the condition 
that different parameters are varied. We then evaluate the 
performance of the system employing different beamformers 
as well as a linear or a ML receiver in an environment in which 
the channels are correlated. We note that for a linear receiver, 
the major computation occurs in the inversion of the Toeplitz 
matrix for which the complexity is of order 0{LM) [61], 
where L is the length of signal vector and M is the number of 
transmitter antennas for the MISO system. On the other hand, 
the complexity using a ML detector for this MISO system 
transmitting the Toeplitz code is of order /i^^ where ji is 
the constellation cardinality. Thus, for a reasonably large con- 
stellation and/or a comparatively large number of transmitter 
antennas, the ML receiver is substantially more complex than 
a linear receiver Finally, we compare the performance of the 
Toeplitz code with some other known efficient codes. 
Example 1: In this example, we examine the performance 
of Toeplitz STBC for a MISO communication system with 
independent channel fading, i.e., S = I. The system is 
equipped with a linear ZF detector at the receiver end. For the 
Toeplitz STBC, we choose B = I in Eq. ( l30l l. The following 
three experiments are performed: 

1 ) We fix the number of transmitter antennas to be i\/ = 4 
and the length of signal vector s to be L = 8, and the 
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Fig. 2. The average BER performance of the proposed ToepHtz STBC for 
different L. 
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Fig. 3. The average BER performance of the proposed Toeplitz STBC for 
the MISO system with different number of transmitter antennas. 



symbol transmission data rate is therefore Rg ~ L/N = 
0.7273 symbols pcu. The signals are randomly selected 
from the different constellations of 4-QAM, 16-QAM and 
64-QAM. The signals are transmitted through the MISO 
system having zero-mean unit-variance i.i.d. Gaussian 
channels and additive white Gaussian noise as described 
in Section |II] and the BER curves are plotted in Fig. [T] in 
which the three different curves correspond to the perfor- 
mance of the system using the three signal constellations 
respectively. It should be noted that different constellation 
size results in different transmission bit rates. For the 
system we examine, the bit rates are i?t = 1.4545, 
2.9091, and 5.8182 bits pcu corresponding to 4-QAM, 
16-QAM and 64-QAM respectively. Therefore, for larger 
constellation size, worse BER performance is expected. 
This is indeed the case as shown by the three BER curves 
plotted in Fig.[T]from which it is observed that, for a BER 
of 10"^, the difference in SNR between 4-QAM and 16- 
QAM is approximately 3dB and that between 16-QAM 
and 64-QAM is approximately 5 dB. 
2) In this experiment, we fix the signal constellation to be 
16-QAM for a four transmitter antenna MISO system 
and perform simulations for different signal lengths, 
L ~ 4, 8, 16, 32. For the different choices of L, the 



system has different transmission symbol rates, which 

are Rs = 0.5714, 0.7273, 0.8421, 0.9143 symbols 

pcu, respectively. The channel and noise assumptions 

are the same as those in the previous experiment. The 

BER curves at different SNR are plotted in Fig. |2l It 

can be observed that the longer is the transmitted signal, 

the worse is the system BER performance. Again, this 

is due to the fact that larger L corresponds to higher 

transmission date rate resulting in worse performance. 

3) In this experiment, we vary the number of the transmitter 

antennas AI. An increase in M increases the diversity 

and decreases the transmission symbol rate. Therefore, it 

is expected that the performance of the system will be 

enhanced with the increase of the number of transmitter 

antennas. This is indeed the case as illustrated in Fig. [3] 

where we compare the performance of the MISO systems 

having M = 2,4, and 8 antennas with L = 8 in the 

Toeplitz codes and the signals selected from a 16-QAM 

constellation. 

Example 2: In this example, we test the performance of 

Toeplitz STBC for correlated channels in a MISO system 

equipped with four transmitter antennas in a linear array 

and the one receiver antenna on the normal to the axis of 

the transmitter antenna array ("broadside"). For small angle 

spread, the correlation coefficient between the mith and m2th 

transmitter antennas is [43], [62] 

1 /■'" 



2tt 



-j27r(mi - TO2)A— sin6' ) dO (79) 



where dt is the antenna spacing, <; is the wavelength of the 

(narrowband) signal, and A is the angle spread. Here in our 

simulations, we choose dt = 0.5"; and A = 5°. We examine 

the performance of the MISO system transmitting 4-QAM 

signals in the following three cases using Toeplitz STBC 

having the structures; 

i) B = -Tfj'^- This is an approximately optimal transmission 

matrix at sufficiently high SNR in the minimization of the 

Chernoff bound under the assumption that B is a square 

matrix. The approximate optimality can shown as follows: 

For Eq. ( |76] |, under high SNR, we ignore the identity 

matrix I in the denominator and obtain 

(8a2)A^ 



P(s^s'|B) < 



2dct(EA'^(e)A'B(e)) 

(8^2 )A^ 



2 det(S)det(B^r«(e)r(e)B)^^^^ 

To minimize the right side of Eq. dSOl l, we maximize the 
second determinant in the denominator. We note that 



det(B^r^(e)T(e)B) 



M 



< 



det(T«(e)T(e))n [BB^] 



m—1 



<det(T«(e)r(e))(^l:tr(BB^; 
dct(r^(e)r(e)) 
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M^-f 



(81) 
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Bop. This minimizes Chernoff bound on PEP as 



The first inequality is due to Hadamard Inequality [44] 
and equality holds iff BB^ is diagonal. The second 
inequality is due to geometric mean being no larger than 
arithmetic mean, with equality holding iff BB^ has equal 
diagonal elements. Finally, the trace of BB^ is equal to 
unity due to the power constraint. Hence, the condition 
for maximum in Eq. dSTT i is that B is a scaled unitary 
matrix of which B = I is one choice, 
ii) B = Bop. This is the optimal solution to minimize the 
worst case PEP derived in Theorem|2]and can be obtained 
numerically by solving the convex optimization problem 
of Eq. ( l64l i and then using the result in Eq 
iii) B 

described in Corollary [T] 
In our simulations here, the transmitted signal vector is of 
length L = 10 and each of the symbols is randomly selected 
from the 4-QAM constellation. At the receiver, the signals in 
each of the three cases are detected separately by a linear ZF 
detector and a ML detector and the respective performances 
of the two detectors are examined. (As mentioned in the end 
of the last section, we cannot obtain an exact optimum B 
for the ZF receiver Nonetheless, we will employ the two 
optimum transmission matrices derived for the ML receiver 
Bop and Bop to the case of ZF receiver to see how the 
performance is improved). We note that due to the variation 
of the channel fading, the dimension K of the optimum 
transmission matrices in Cases ii) and iii) change with SNR. 
For the specific correlated channel described in Eq. ( |79] l, for 
Case ii), we found that K = 1 when SNR < 8dB and K = 2 
at higher SNR, whereas for Case iii), we found that K = 1 
when SNR < lOdB and K = 2 at higher SNR. Therefore, the 
transmission data rate for these two cases is Ri = 
symbols pcu, which is lower than i?2 = ^ ,'^^_-^ 
of B = Ij\f . For a fair comparison, we choose two different 
structures for B in Case i) in the following experiments; 
1) We maintain the transmission data rate in Case i) the 
same as that in case iii). This is realized by setting 
B = [ -^1k,0{k,m-k)\ in Case i) with K being 
the dimension of B^Bop. We evaluated the error per- 
formances of the systems equipped with different B in 
all three cases and the results are shown in Fig. |4] from 
which the following observations can be made: 

• For the system employing a ML detector, perfor- 
mance of Case ii) and iii) are superior to that of Case 
i), confirming the theoretical analyses in Theorem |2] 
and Corollary [T] 

• For the system employing a ML detector, the BER 
performance for Cases ii) and iii) employing Bop 
and Bop respectively are very close. This shows 
that Chernoff bound is tight for this system. Close 
performance in the two cases is also true for the 
system using a ZF detector. 

• Although Bop and Bop are optimal transmission ma- 
trices developed for the ML detector, they are equally 
effective in providing substantial performance im- 
provement for the same system employing a linear 
ZF detector 



K 
L+K-1 

for the case 
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Fig. 4. The average bit error rate comparison of the proposed Toeplitz STBC 
with i) B = [Ik, Oa/— k], ii) Bop and iii) Bop. The performances are shown 
for both linear ZF detectors and ML detectors. 



!\/l=4 



10 



H 10" 

DC 



.■= 10" 



10" 



10" 





- * - ZF, B=l 
A ZF, case iii) 
□ iUlU B=i 

^ ° X 

AD « 
N. AD * 

X AD 

\ An 



10 15 

SNR [dB] 



20 



25 



Fig. 5. The average bit error rate comparison of the proposed Toeplitz 
STBC with B = Ijv/ and Bop. The performances are shown for both linear 
ZF detectors and ML detectors. 

• At lower SNR, we have K ^ 1, i.e., only one 
transmitter antenna is effective. Therefore, given a 
coded system, Hnear ZF and ML detectors provide 
the same performance. 

2) In the second part of the experiment, we put B = Ij\/ 
in Case i) and examine its performance. The error per- 
formance for such a choice is shown in Fig. |5] Here, 
the system using B = Im has higher transmission data 
rate than those in Cases ii) and iii). For the sake of 
comparison, we have re-plotted in Fig. |5] the performance 
curves from Fig. |4]of Case iii) corresponding to the uses 
of Bop as a transmission matrix. (Since the performance 
of Cases ii) is almost the same as that of Case iii), we 
have omitted here the performance curves corresponding 
to the use of Bop). It should be noted that when the 
signals are detected by a ML detector, the system coded 
with B = Im has higher diversity gain over the system 
with Bop. This is due to the fact that water-filling strategy 
may not employ all the available transmitter antennas 
for correlated channels. Specifically in this example, the 
effective number of antennas for Bop is K < 2 < 4. 



However, the optimal coding gain achieved by Bop 



with 



ML detectors ensures a better performance. It is also im- 
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Fig. 6. The average bit error rate comparison of the proposed Toeplitz STBC p^g 7 ^he average bit error rate comparison of the proposed Toeplitz STBC 



with other STBC of unit rate. 



with the orthogonal STBC. 



portant to note that the employments of B = I]^f and Bop 
result in a relatively large difference in performances, 
revealing that the upper bound on PEP given in Eq. dSOl l is 
not tight. Thus, even though this bound is quite commonly 
employed in STBC designs for independent channels, the 
results here show that this relaxed bound is a poor design 
criterion for an environment of highly correlated channel 
coefficients. 
Example 3: In this example, we compare the BER perfor- 
mance of Toeplitz STBC with other STBC for independent 
MISO channels. Here again, we choose B = I for Toeplitz 
STBC. The experiments are performed for the two cases in 
which the number of transmitter antennas in the communica- 
tion system are M = A and 7\f = 8 respectively: 

1) Af = 4 transmitter antennas and a single receiver antenna: 
We compare BER performance of Toeplitz STBC with 
other rate one STBC [29], [34], [36], [63]: 

• Quasi-orthogonal STBC. The code for four transmit- 
ter antennas was presented in [34], and the maxi- 
mization of its coding gain was subsequently shown 
in [36]. 

. Dense full-diversity STBC [29] 

. Multi-group decodable STBC [63] 
For the Toeplitz STBC, we choose L = 9 for which 
the symbol transmission data rate is Rg = L/N = 3/4 
symbols pcu. To achieve a fair comparison, the same 
transmission bit rate is imposed on all the codes such 
that signals are selected from 256-QAM constellation for 
ToepHtz STBC and from 64-QAM for the other full-rate 
STBC. Therefore, the same transmission bit rate, Ri, ~ 6 
bits pcu, is employed for all the systems. At the receiver, 
the Toeplitz STBC is processed by a linear ZF equalizer 
followed by a symbol-by-symbol detector For the other 
full-rate STBC, we examine the two cases in which the 
signals are processed by a) a ML detector and b) a 
linear ZF receiver The BER curves are plotted in Fig. |6] 
When a linear ZF equalizer and a symbol-by-symbol 
detector is applied at the receiver, it can be observed that 
Toeplitz STBC outperforms "quasi-orthogonal" STBC 
and "dense" STBC, and at higher SNR, its performance 



is superior to multi-group code. It is also interesting to 
observe that at higher SNR, for Toeplitz STBC with linear 
ZF receivers, the performance is also superior to that of 
the Multi-group STBC using a ML receiver. In fact, for 
the range of SNR tested, the slope of its BER curve is 
the same as those of the "dense" STBC and the "quasi- 
orthogonal" STBC processed by ML detectors, indicating 
they have the same diversity gain. 
2) We now consider the system having M ~ 8 transmitter 
antennas. For the Toeplitz code, we choose L = 35 and 
therefore, the symbol transmission data rate is L/N = 
5/6 symbols pcu. We compare the bit error rate perfor- 
mance of our Toeplitz code with that of the orthogonal 
STBC having symbol transmission rate of: 
i) 1/2 symbols pcu [2], [3], [7] and 
ii) 5/8 symbols pcu [64] (this the highest symbol rate 
achievable by the orthogonal STBC applied to an eight 
transmitter antenna system). 

To achieve a fair comparison, the transmitted signals are 
selected from a 64-QAM constellation for our Toeplitz 
code, a 256-QAM constellation for the 5/8 rate orthogo- 
nal code and a 1024-QAM constellation for the 1/2 rate 
orthogonal code. Hence, all of the codes have the same 
transmission data rate in bits, i.e., R^ ~ 5 bits pcu. At the 
receiver end, the orthogonal STBC is decoded by a linear 
ZF detector for which, because of the orthogonality, the 
performance is the same as that of a ML detector For 
Toeplitz STBC, the signals are decoded separately by a 
linear ZF receiver and a ZF-DFE receiver. The average 
bit error rate for these codes are plotted Fig. [T] It can 
be observed that the performance of the Toeplitz code 
detected with a linear ZF receiver is superior to that of 
the |-rate orthogonal STBC when the SNR is less than 
or equal to 25 dB. When the Toeplitz STBC is received 
by a ZF-DFE receiver, due to the higher coding gain, 
its performance is significantly better than that of the 
orthogonal STBC. In Fig. |7] at 10~^ the Toeplitz code 
with a ZF-DFE receiver outperforms the orthogonal code 
by about 4 dB. 
It should be noted that for the Toeplitz code, both lin- 
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Fig. 8. The average symbol error rate comparison of the proposed Toeplitz 
STBC with the orthogonal STBC. 



ear ZF and ZF-DFE receivers can achieve full-diversity. 
However, from Fig. |7] while the slope of BER curve 
for Toeplitz code with ZF-DFE receiver is similar to 
those of the orthogonal codes, the slope of the curve 
for the Toeplitz code with linear ZF receiver is not as 
steep. Recall that the diversity gain of a communication 
system is defined at high SNR and here, the upper end 
of the SNR range is not sufficiently high. To show full 
diversity for both systems, we need BER at higher SNR, 
the evaluation of which demands exorbitant computation 
for the parameters in this example. To circumvent this 
difficulty, we choose to compare the symbol error rate 
(SER) obtained by the use of the Toeplitz code with 
linear ZF receiver to that obtained by the use of the 5/8 
orthogonal code. The results are shown in Fig. [8] from 
which it can be observed that the two SER curves have 
the same slope for SNR above 30dB, indicating the same 
diversity gain for both codes. Thus, we can see that the 
Toeplitz code with a linear ZF (or more sophisticated) 
receiver indeed achieves full diversity. 

VIII. Conclusion 

In this paper, we have presented a general design criterion 
for full-diversity linear STBC when the signals are transmitted 
through a MISO communication system and processed by 
a linear receiver This is, to our knowledge, the first de- 
sign criterion for linear receivers to achieve full diversity. 
Specifically, we proposed a linear Toeplitz STBC for a MISO 
channel which satisfies the criterion and achieves full-diversity. 
We have shown that such a code possesses many interesting 
properties, two of which recapitulated here are of practical 
importance; 

1) The symbol transmission rate for the code approaches 
one when the number of channel uses {N > M) is large. 

2) If the signalling scheme has a constellation for which the 
distance between the nearest neighbours is nonzero (such 
as QAM), then employing the Toeplitz code results in a 
non-vanishing determinant. 

When employed in a MISO system equipped with a linear 
receiver (ZF or MMSE), the Toeplitz code can provide full 



diversity. Furthermore, when the number of channel uses is 
large, in an independent MISO flat fading environment, the 
Toeplitz code can approach the Zheng-Tse optimal diversity- 
multiplexing tradeoff. 

When employed in a MISO system equipped with a ML 
detector, for both independent and correlated channel coeffi- 
cients, we can design the transmission matrix inherent in the 
proposed Toeplitz STBC to minimize the exact worst case 
average pair-wise error probability resulting in full diversity 
and optimal coding gain being achieved. In particular, when 
the design criterion of the worst case average pair-wise error 
probability is approximated by the Chernoff bound, we obtain 
a closed-form optimal solution. 

The use of the Toeplitz STBC (having an identity trans- 
mission matrix) in a MISO system fitted with a ZF receiver 
has been shown by simulations to have the same slope of 
the BER curves to other full rate STBC employing a ML 
detector, whereas even better performance can be achieved 
by using receivers (such as ZF-DFE) more sophisticated than 
the linear ones to detect the Toeplitz code. For correlated 
channels, employing the optimum transmission matrices in the 
Toeplitz code results in substantial additional improvements 
in performance to using the identity transmission matrix. This 
substantial improvement of performance is observed in either 
case for which a ML or a ZF receiver is used. 
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